Dense Packings of Congruent Circles 
in Rectangles with a Variable Aspect Ratio 



Boris D. Lubaclievsky Ronald Graham 

bdl@bell-labs. com graham@ucsd. edu 

Bell Laboratories University of California 

600 Mountain Avenue at San Diego 

Murray Hill, New Jersey La JoUa, California 



Abstract 

We use computational experiments to find the rectangles of minimum area into 
which a given number n of non-overlapping congruent circles can be packed. No as- 
sumption is made on the shape of the rectangles. Most of the packings found have 
the usual regular square or hexagonal pattern. However, for 1495 values of n in the 
tested range n < 5000, specifically, for n = 49, 61, 79, 97, 107, ...4999, we prove that 
the optimum cannot possibly be achieved by such regular arrangements. The evidence 
suggests that the limiting height-to-width ratio of rectangles containing an optimal 
hexagonal packing of circles tends to 2 — \/3 as n — > oo, if the limit exists. 

Key words: disk packings, rectangle, container design, hexagonal, square grid 
AMS subject classification: primary 52C15, secondary 05B40, 90C59 

1. Introduction 

Consider the task of finding the smallest area rectangular region that encloses a given 
number n of circular disks of equal diameter. The circles must not overlap with each other 
or extend outside the rectangle. The aspect ratio of the rectangle, i.e., the ratio of its height 
to width, is variable and subject to the area-minimizing choice as well as the positions of 
the circles inside the rectangle. 

Packing circles in a square has been the subject of many investigations [(jLj . |N()lj . 
NU2j . |NU3j . Because the aspect ratio is not fixed in our present problem, the solutions 



are typically different from the dense packings in a square. For example, the density 7r/4 
0.785... of the proved optimum (see |N03j) for packing 25 congruent circles in a square (see 
Figure ITTk ) can be increased to 257r/ (26(2 + v^)) = 0.809..., if we let the rectangle assume 
its best aspect ratio (see Figure ITTTb ). Our experiments offer only three values of n for which 
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Figure 1.1: The best packings found for 25 circles a) in a square, and, b) in a rectangle with 
variable aspect ratio. 




Figure 1.2: A 224-circle fragment of the packing with the highest density found in a long 
rectangle with a fixed aspect ratio. 
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the dense packing in a square is also a solution to our present problem: n = 4, 9, and of 
course, n = 1. 

Similarly, dense packings in long rectangles with a fixed aspect ratio usually do not yield 
solutions for our problem. According to a long-standing conjecture attributed to Molnar 
(see IjFiiredij ). dense packings in long fixed rectangles tend to form periodic up and down 
alternating triangular "teeth," as in the example shown in Figure 11.21 We can usually 
increase the density of such packings by slightly changing the aspect ratio of the rectangle. 

The present problem has been occasionally mentioned among packing problems ; for 
example, Web enthusiasts have recently began discussing it. However, the problem was 
considered as early as in 1970 in jEHSi, where the optimum packings were determined 
for all n < 8 and conjectured for 9 < n < 12. We learned about Ruda's results after 
we performed our computations and made our conjectures. Fortunately, our conjectures 
restricted to the values n < 12 agree with the results and conjectures in |R,udaj ! 

It appears that the complexity of the proofs of optimality in the packing-circles-in-a- 
square problem increases exponentially with n; computer generated/assisted proofs of opti- 
mahty of such packings do not reach n = 30 (see |NU3j ; proofs that do not utilize computers 
have not been found except for quite small values of n), while conjectures extend to n = 50 
and even much larger values (see |NDRJ). 

Similarly, it seems to be difficult for the present problem to prove optimality for the 
configurations we find. This paper describes an experimental approach, where good packings 
are obtained with the help of a computer. Some of these packings hopefully will be proved 
optimal in the future. At present, most of the statements about the packings made in this 
paper are only conjectures, except for a few which we explicitly claim as proven. For instance, 
we prove that the best packings found are better than any other in their class, e.g., that 
the packing of 79 circles with a monovacancy (see Figure , while non-optimal, is better 
than any hexagonal or square grid packing of 79 circles without a monovacancy. Also, we 
prove that n = 11 is the smallest n for which a hexagonal packing as in Figure ITTl is better 
than any square grid packing of n circles. 

2. A priori expectations and questions about best packings 

Since it is well known that the hexagonal arrangement of congruent circles in the infinite 
plane has the highest possible density, (see |Thj . |Ft] . [FG],[01erj), before plunging into our 
experiments, we expected to obtain good finite packings by "carving" rectangular subsets 
out of this infinite packing. It will be useful to classify here such finite arrangements. For 
that we will discuss packings in Figures ITTTb . 12. H 12.2b . HT^h . and 12.3b . irrespective of their 
optimahty (which will be discussed in the following sections) . For the purpose of exposition, 
we will pretend in this section that there are no holes in the arrangements in Figures f2.2b 
and 12.3b . that is, the question mark in the figures is covered with an additional circle of the 
common radius. 

Under such an assumption, the arrangement in Figure 1?!^ consists of h = 5 alternating 
rows with w = 16 circles in each, a total of n = w x h = 80 circles. Arrangements as in 
Figures ITTTb . 12. 1112. 2b . and 12.2b . are obtained by depleting such w x h arrays of some circles 
along one side of the array. The ways of performing the depletion depend on the parity of h. 
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Figure 2.1: The best packing found for 11 circles in a rectangle. 




Figure 2.2: Packings of 49 circles in a rectangle: a) the best in the class of hexagonal 
packings, b) a best in the class of hexagonal packings with monovacancies (one of 17 equally 
dense packings with the hole), c) the best we found. 
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Figure 2.3: Packings of 79 circles in a rectangle: a) a best in the class of hexagonal packing 
with possible monovacancies, and, b) the best we found. 




If h is even, then the depletion can be done in only one way by removing h/2 circles. 
Figure ITTTb presents an example with h = 2. 

The case of an odd h presents two possibilities: we can remove [_h/2\ circles on a side as 
in the example of Figure or \_h/2\ +1 circles on a side as in the example of Figure EI2t- 

Alternatively, one can consider rectangular square grid arrangements as candidates for 
the best packings, e.g., see those in Figure |2]U The density of a square grid arrangement 
is fixed at 7r/4 independent of n. When both sides h and w of the carved rectangle tend to 
infinity, the hexagonal packing density tends to ^T/(2^/3) which is the density of the hexagonal 
packing in the infinite plane. Since 7r/(2v^) > 7r/4, one naturally wonders for which n the 
best hexagonal grid arrangement becomes better than the square grid arrangement. 

A natural question is whether or not both arrangements exhaust all possibilities for the 
best packing. In other words, does there exist an optimal packing of n disks in a rectangle 
such that it cannot be represented as being carved out by a rectangle from either square grid 
or hexagonal grid packing of the infinite plane? 

3. Results of compactor simulations 

To tackle the problem by computer, we developed a "compactor" simulation algorithm. 
The simulation begins by starting with a random initial configuration with n circles lying 
inside a (large) rectangle without circle-circle overlaps. The starting configuration is feasible 
but is usually rather sparse. Then the computer imitates a "compactor" with each side of 
the rectangle pressing against the circles, so that the circles are being forced towards each 
other until they "jam." Possible circle-circle or circle-boundary conflicts are resolved using 
a simulation of a hard collision so that no overlaps occur during the process. 

The simulation for a particular n is repeated many times, with different starting circle 
configurations. If the final density in a run is larger than the record achieved thus far, it 
replaces this record. Eventually in this process, the record stops improving up to the level of 
accuracy induced by the double precision accuracy of the computer. The resulting packing 
now becomes a candidate for the optimal packing for this value of n. 

In this way we found that the square grid pattern with density 7r/4 supplies the optimum 
for n = 1, 10, and that ?7, = 11 is the smallest number of circles for which a density better 
than 7r/4 can be reached; the density of this packing is ll7r/(16(l + v^) = 0.790558.. The 
corresponding conjectured optimum pattern for n = 11 is shown in Figure lTTl The pattern is 
hexagonal. For n = 12 and 13 the square grid pattern briefly takes over as the experimental 
optimum again (see Figure 1231 for the case of n = 12). Remarkably, Ruda |Rudaj proves 
that for n < 8, the square-grid packings are optimal. His conjectures for 9 < n < 12 also 
agree with our findings. 

Our simulation results show that for n > 14, the best densities are larger than 7r/4. This 
statement is also easy to prove, considering examples of (not necessarily optimal) two-row 
hexagonal packings like that in Figure II. lb for odd n = 2m + 1 > 14, or the ones with 
equal length rows for even n = 2m > 14. Assuming the circles have radius 1, the following 
inequalities 




(1) 
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Figure 3.1: Two equally dense best packings found for 15 circles in a rectangle. 



for n = 2m + 1 and 

(2m+ 1)(2 + V3) < 8m (2) 

for n = 2m have to be satisfied. The left-hand sides in and (j21) are the areas of the 
enclosing rectangles for the two-row hexagonal packings, and the right-hand sides are the 
areas for the corresponding square grid packings. It is easily seen that all integers m > 7 
satisfy either inequality. These correspond to all integer values of n > 14. 

A surprise awaited us at the value n = 15. We found two different equally dense packings. 
One, shown in Figure ITTk . is of the expected hexagonal type. However, the other, shown in 
Figure ITTb . is not, nor can it be carved out of the square grid. It is easy to verify that two 
equally dense packings a and b of these types also exist for any n of the form n = 15 -|- Ak, 
k = 1,2,3,.... Packing a, such as that in Figure IXTk . has two alternating rows, with one 
row one circle shorter than the other, and with the longer row consisting of w = 8 + 2k 
circles. Packing b, such as that in Figure ITTb . has 4 rows, the longest row having w = A + k 
circles, three bottom rows alternate, the middle one having w — 1 circles, and the 4th row 
of w circles stacked straight on the top (or on the bottom; these would produce the same 
configuration). Our experiments suggest that for k = 1 and 4, i.e., for n = 19 and 31, these 
pairs of configurations might also be optimal; however, they are provably not optimal for 
other values of > 0. The packings of pattern b for n = 15, 19, and 31, if they are proved 
to be optimal, answer positively the second question posed in Section |21 

Thus, the optimal container for 15 bottles apparently can have two different shapes, 
as can the optimal containers for 19 and 31 bottles! In an obvious way, more than one 
optimum container shape also exists for square grid packings of n circles if n is not a prime. 
See Figure 12.41 for the example of = 12; there are three different rectangles which are 
equally good and probably optimal for n = 12. Apparently, three is the maximum number 
of rectangular shapes, any number n of circles can optimally fit in. It seems that for n = 
4, 6, 8, 9, 10, 15, 19 and 31, there are exactly two shapes. For all other n tested, only one best 
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aspect ratio was found experimentally. 

Another packing surprise awaited us at n = 49. Figures f2.2fe i and 12.2b show two among 
several equivalent packings achieved in our simulation experiments. The configuration a is 
hexagonal, while configuration h contains a monovacancy, that is, a hole than can accom- 
modate exactly one circle. Monovacancies in hexagonal arrays of congruent circles appear 
often in the simulation experiments for 14 < n < 49 but only in packings of inferior quality, 
so that a better quality hexagonal packing without monovacancies could always be found. 
No higher density hexagonal packing without monovacancies was found for = 49. 

Such large n already present substantial difficulties for our simulation procedure. The 
procedure failed to produce a packing which is better than those in Figures 12.2b and 12.2b . 
However, it is easy to prove that the density of a hexagonal packing in a rectangle with 
a monovacancy can be increased. For example, to improve the packing in Figure 12.2b we 
relocate circle B into the vacancy and then rearrange circles A and C along the side. This 
reduces the width of the rectangle by a small but positive 5 as shown in Figure 1?!^ . where 
5 = 2- \f%Jl = 0.13879.. of the circle radius. The density of c is 497r/(2(l + V^){3A -6)) = 
0.83200266... 

It is not known whether or not the resulting configuration c can be further improved. 
However, using the exhaustive search we show (see next section) that the optimum packing 
for 49 circles in a rectangle, whatever it might be, cannot be purely hexagonal. The value 
n = 49 is apparently the smallest one for which neither square-grid nor perfect hexagonal 
pattern delivers the optimum. 

4. Exhaustive search 

The simulation method becomes progressively slower for increasing n. However, by 
exercising the simulation for smaller n we are able to refine the idea of the class of packings 
which might deliver the optimum. We have chosen the class that consists of square grid and 
hexagonal packings and their hybrids; we also allow for monovacancies in configurations of 
the class, although these configurations are never optimal. We have extended our computing 
experiments to larger values of n using exhaustive search. For a given n, the method simply 
evaluates each candidate in the class by computing the area of the rectangle and selects the 
one with the smallest area. Note that for each n in this class, there are only a finite number 
of candidates. 

A general packing in the class (see Figure ^HJ, consists of h + s rows and has d mono- 
vacancies. The h rows are hexagonally alternating, and the s rows are stacked directly on 
top of the previous row as in square grid packings. The longest row consists of w circles. 
Assuming all monovacancies are filled with circles, among the h hexagonal rows, h- rows 
consist of w — 1 circles each, the remaining h — h^ rows consist of w circles each, and among 
the s square grid rows, s_ rows consist of w — 1 circles each, and the remaining s — s_ rows 
consist of w circles each. The number of circles in the configuration is 



Table 14.11 lists the best packings found of n circles in a rectangle with variable aspect 
ratio for n in the range 1 <n < 53, and Table W7I\ continues this list for 54 < < 213. The 




(3) 
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Figure 4.1: A "general case" member of the class of packings among which we search for the 
optimum; here w = 5, /i = 5, /i_ = 2, s = 2, s_ = 1, and d = 3. 

packings in Table 14.11 are obtained by the simulation procedure described in Section El and 
verified by the exhaustive search. Most packings in Table 14.21 were generated only by the 
exhaustive search. An entry in either table consists of the n, followed by the set of integers 
which represent the parameters of the packing structure as explained above: 
w, the number of circles in the longest row, 

h, the number of rows arranged in a hexagonal alternating pattern, 
h-, the number of rows that consist of w — 1 circles each, 

s, the number of rows, in addition to h rows, that are stacked in the square grid pattern. 
Column s is absent in Table 14.21 because, as explained in Section El no optimum packing 
found for n > 31 has s positive. Also, for any n, no optimum packing was found with s_ > 0; 
thus, the S- column is omitted. 

In most packings presented in these two tables, the number of monovacancies d = 0, and 
so the d column is omitted. A few entries where monovacancies are possible are marked with 
stars. The non-starred entries describe the best packings found. The entries marked with 
stars cannot be optimal, as explained in Sectional 

An improved packing for the marked case of 49 circles can be obtained as described in 
Sectional (see Figure . The next marked case n = 61 is similar; Table l^?T1 lists the variant 
a for it with h = 3 and h_ = \_h/2\ +1 = 2. The following marked case is n = 79. Here 
too we improve the packing with the hole by relocating side circles. A, B, C, D, and E as 
shown in Figure EISl where the value of improvement S is the same as in the case n = 49. 

This method of relocation applies to all marked cases in the tables when h is odd. For 
example, for the case of h = 5, when /i_ = 3 the relocation can be done as illustrated 
in Figure 14.31 This applies to the cases n = 97, 107, and 142 in Table 14.21 In all these 
cases, before using this method, we should replace listed in the tables variant that has = 
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Table 4.1: Packings of n circles in a rectangle for 1 < n < 53 with w circles in a row, h 
hexagonal rows, h- of which consist of w — 1 circles each, and s square grid rows of w circles 
each. Except for the case n — A9 marked with a star, all packings are the best we could find. 
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\h/2\ = 2 with the variant that has /i_ = [_h/2\ +1 = 3. This would produce a configuration 
with the pattern of the side as in Figure 14.3b . We would then improve this configuration by 
relocating circles A, B, C, D, and E to yield a configuration as in Figure HTHb . The value 

of the improvement here is <5 = 2 - 0.5v^ - 3^/^(273 - l)/{2\/4: - Vs) = 0.05728... of the 
circle radius. A similar method works for even h. 

Sometimes during these improvements, some circles become the so-called rattlers, i.e., 
they become free to move and hit their neighbors, like the two unshaded circles in Figure EHl 
For glass bottles tightly packed in an empty box, packings with rattlers should definitely be 
avoided even though the box area was minimal! 

The case n = 79 and several others are marked in Table 14.21 with two stars to signify 
that the packing not only may contain a monovacancy, but in fact, it must: any hexagonal 
packing of n circles without a monovacancy is provably worse than the one represented in 
the table when n is marked with two stars. 

All such cases in the table also happen to have h_ = 0. Also, a double-star marked 
entry in Table 14.21 always has a discrepancy between the n computed from the parameters 
of the packing structure using (jni),(here it would be n = wh), and actual n. In all the other 
cases, the n computed from the packing structure parameters using Q will always match the 
correct n because there are no monovacancies. However, no entry without a monovacancy is 
possible in the double-star marked cases. Another observation: an entry n + 1 that follows 
an entry n marked by one or two stars, always has the same w and h, and hence its packing 
fits in the same rectangle as that of the entry n. 

5. Best packings found for larger n 

The exhaustive search procedure of Section |3] produces the best packings in the class 
defined in Section |3 for values of n on the order of several thousands. The features we 
observed for n < 213 hold until n = 317. Namely, only one best size rectangle exists for each 
n. Most of the best packings in the search class are perfectly hexagonal and unique for their 
n. Each such packing can be described by the set w, h, and with the latter parameter 
taking on up to three possible values: 



A few exceptional cases have a single monovacancy. These are similar to the one-star marked 
cases with odd h and with two options for (the choice of the option defines the presence 
or absence of the vacancy), or they are similar to the two-star marked cases with = 0. 

The case n = 317 deviates from this pattern. Here, the best packing in the class has 
parameters w = 27, h = 12, h_ = 6, and d = 1. Unlike the one-star marked cases, h is even, 
and unlike the two-star marked cases, /i_ is non-zero. Such a new type of packing recurs 
with one hole for n = 334 [w = 34, h = 10, and h_ = 5) and then for n = 393 {w = 40, 
h = 10, and /i_ =5), but in the latter case with two monovacancies. In other words, the 
best possible packing for 393 circles in a rectangle in the class of hexagonal packings with 
possible monovacancies must have two of them. 
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15 


9 


4 


171 


16 


11 


5 


*211 


24 


9 


5 


92 


23 


4 





132 


22 


6 





172 


25 


7 


3 


212 


24 


9 


4 


93 


19 


5 


2 


133 


27 


5 


2 


173 


35 


5 


2 


213 


36 


6 


3 



Table 4.2: Packings of n circles in a rectangle for 54 < n < 213 with w circles in a row, h 
hexagonal rows, h_ of which consist of w — 1 circles each. Except for the cases marked with 
stars, all packings are the best we could find. 
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Figure 4.2: Numbers of circles n{i) for which the optimum packing in a rectangle is not 
perfectly hexagonal for 49 < n{i) < 5000 plotted versus the rank i. 
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Figure 4.3: Improving a packing with a monovacancy when h — 5, — 2: a) an original 
hexagonal packing with a monovacancy, b) the transformed packing. 
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With two holes, we have more freedom to improve the quahty of the packing than having 
just one hole. Several ways are possible, including combining the ones described before for 
cases of smaller n. Same as for the other packings with monovacancies, the optimal packing 
for n = 393 is not known, but the exhaustive search in our experiments proves that it cannot 
be a hexagonal packing. 

Naturally, the best packing (in the search class) of n = 394 circles also has parameters 
w = 40, /i = 10, /i_ = 5 and d = 1, i.e., a single monovacancy. The same rectangle also 
accommodates the best packing in the class of n = 395 circles, with the same w = 40, h = 10, 
h_ = 5. Since the latter packing is without holes, it is conceivable that it is optimal. 

The next outstanding case is n = 411 where w = 38, h = 11, h_ = 6, and d = 1. This 
is the smallest n where the h is odd and h- = [h/2\ +1, that is, the packing looks like the 
one in Figure I2I2K, but unlike the latter it has a hole. An equivalent packing that looks like 
the one in Figure HT^h . where = \_h/2\ = 5 also exists, but unlike the latter it has not 
one but two holes. For next value n = 412, we have a standard one-star marked situation 
with the same sizes of the rectangle, i.e., w = 38, h = 11 and /i_ = [/i/2j = 5 for the variant 
with one hole and h- = lh/2\ +1 = 6 for the variant without holes. 

The smallest n for which as many as three monovacancies exist in the best packing in 
the class is n = 717. The corresponding packing parameters are w = 48, h = 15, = and 
no hexagonal packing with a smaller number of holes can be better or even as good as this 
packing. Similarly, with four monovacancies, the smallest n is n = 2732 [w = 86, h = 32, 
and h_ = 16), and for five monovacancies, n = 2776 is the smallest {w = 103, h = 27, and 
h- = 0). No optimum packing in the class for n < 5000 has six or more monovacancies. 

In all, we found 1495 values of n on the interval 1 < n < 5000 for which the best packing 
in the class must or can have monovacancies. As explained above, for each such n the best 
packing provably cannot be of a pure square-grid or hexagonal pattern. It is not proven 
though, but we believe it to be also true, that those 1495 values of n are all such irregular 
values among the considered 5000 values. 

The chance to encounter such an irregular n seems to increase with n. For example, here 
are all the experimentally found irregular values of n among 100 consecutive values on the 
intervals 401 + lOOOfc < n < 500 + lOOOfc for A; = 0, 1, 2, 3, 4: 

17 values for 401 < n < 500: 
409, 411, 412, 421, 422, 433, 439, 453, 454, 461, 463, 467, 471, 478, 487, 489, 499. 

24 values for 1401 <n< 1500: 
1401, 1402, 1405, 1409, 1412, 1414, 1423, 1427, 1429, 1434, 1446, 1447, 1451, 1453, 1457, 
1459, 1466, 1468, 1477, 1483, 1486, 1487, 1489, 1497. 

33 values for 2401 <n< 2500: 
2401, 2402, 2406, 2411, 2419, 2421, 2423, 2428, 2429, 2435, 2437, 2439, 2441, 2443, 2446, 
2452, 2454, 2455, 2456, 2458, 2462, 2467, 2469, 2474, 2476, 2477, 2479, 2481, 2487, 2491, 
2493, 2495, 2497. 

33 values for 3401 <n< 3500: 
3407, 3409, 3411, 3412, 3414, 3415, 3418, 3421, 3425, 3428, 3431, 3433, 3436, 3442, 3446 
3447, 3453, 3455, 3459, 3461, 3464, 3467, 3469, 3473, 3476, 3479, 3481, 3487, 3489, 3490, 
3493, 3494, 3499. 

38 values for 4401 <n< 4500: 
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4401, 4404, 4405, 4409, 4411, 4414, 4417, 4419, 4421, 4426, 4430, 4434, 4436, 4438, 4441, 
4443, 4447, 4450, 4453, 4456, 4457, 4458, 4461, 4462, 4467, 4468, 4474, 4476, 4479, 4483, 
4486, 4487, 4491, 4492, 4493, 4495, 4497, 4499. 

Figure 14.21 represents these 1495 irregular values n{i), 1 < n{i) < 5000, beginning 
with n(l) = 49 and ending with n(1495) = 4999, by plotting points {abscissa = i, 
ordinate = n{i)}. 

6. The optimum aspect ratio for large n 

Figure 16.11 contains for each n, a data point with coordinates 

(abscissa = n, ordinate = best aspect ratio found for this n). 

All n < 5000 are represented with the exception of a few n where the best packings found are 
of the square grid type. Also, the hybrid packings like those in Figure ITTb are excluded. In 
this way we assure the uniqueness of the best aspect ratio for each n represented. The aspect 
ratios for the best packings in the search class are described in Section 0] The changes to the 
aspect ratios are due to the S shrinkages of the rectangles in those cases with monovacancies, 
such as those in Figures 1?!^ 12.31 and 14.31 Presumably these are small in number and should 
not be noticeable in Figure lU^Tl 

The data points tend to form patterns of descending and ascending "threads." To ex- 
amine the threads in detail, a rectangular box which is close to the y-axis in Figure 16.11 is 
magnified and represented in Figure 16.21 All data points present in the box can be found in 
Tables ED and 1121 

The steep downward threads correspond to configurations with fixed height, and with 
widths increasing by 1 from one point to the next as we move down the thread. Each such 
thread eventually terminates in either direction which means that the optimum rectangle 
cannot be too fiat or too tall. The less steep upward and downward threads (they are dotted) 
correspond to sequences {w, h), where from one point to the next w increases by 3 or 6 and 
h increases by 1 or 2, respectively. 

To compute the optimum aspect ratio of the rectangle that encloses the optimum packing 
for large n, we analyze the area wasted along the rectangle sides. Note that in the infinite 
hexagonal packing, the uncovered area is s = -\/3 — -k/2 per each 7r/2 of the covered area. 
This is obvious from examining a single triangle XYZ in Figure 16.31 and observing that the 
entire infinite packing is composed of such triangles. The waste s here is the area of the 
central triangle formed by three circular arcs and it is equal to the full area -\/3 of triangle 
XYZ (assuming circle radius is 1) minus tt/2, which is the total area of the three 7r/6 sectors 
covered. (By the way, we remind the reader that the density of the infinite hexagonal packing 
is (v/3 - s)/{n/2) = 7r/{2V^) = 0.90689968....) 

The structure of the uncovered area in the infinite hexagonal packing can be also under- 
stood if we think of each covered circle bringing with it two curved uncovered triangles s, 
those adjacent to this circle on its right-hand side. With such bookkeeping, each triangle 
s will be counted exactly once. We conjecture that these two triangles s per each circle is 
the unavoidable, i.e., fixed, waste in any finite hexagonal packing carved out by a rectangle. 



16 



24 



24. 



27- 



2.9- 



30-' 



32' 



42 



■■■■■44 



52-^^' 



54- 



130 



■68... 



70 



96 



9-9... 



7-3. 



105- 



.7.5. 



..78 

.79-80 

...8-3-^^'" 
8 -5 - ■ ■ ■ ■ ■ ■ 
88 
90 
93 



ll-l 

....1.14- 

....1.17 

..1.20-^^ 



■1-3-7.. 



■■1-80 



■144. 



184 
1-88 



126 

129 
132 
135 



■14-7 1-91192 

1.5.1 1^96 

1-54 

1.5.7158 

141 

165 

172 
175 
179 

181182 



Figure 6.2: A box selected in Figure IFTD magnified. Each data point in Figure IFTTl is replaced 
here ■with the corresponding number of circles. 
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Figure 6.3: Calculation of the waste area adjacent to the sides of the rectangle. 



There will be additional, i.e., variable, waste along the rectangle sides and we are trying to 
minimize this variable waste. 

Along the bottom and top sides, for each 2 units of length, like the side AD of rectangle 
ABCD in Figure the additional waste is the area of rectangle ABCD minus two covered 
quarter-circle areas and minus s. This s represents one of the two curved triangles attached 
to the right-hand side of the circle with the center at B and as an unavoidable waste should 
not be included. Hence the waste per unit length of the top and bottom sides is a = 
(2 - 7r/2 - s)/2 = (2 - V3)/2. 

The waste along the left- or right-hand sides per two alternating rows, i.e., per 2^3 
units of length, is the area of the semicircle that is cut in half by the side plus or minus 
additional area depending on the side. Along the left-hand side, all additional uncovered 
area is additional waste, since we assume that the curved triangles s are attached to the 
right-hand side of the covered circles. The addition consists of two triangles s and two 
halves of triangles s/2 as shown in Figure The left-hand side waste per 2-\/3 of length 
is thus 7r/2 + 3s. Along the right-hand side we subtract from the area of the semicircle the 
outstanding two half triangles s/2 because these are necessary in the infinite packing. The 
right-hand side waste per 2v^ units of length is thus 7r/2 — s. Averaging this for both sides 
and dividing by 2v^, we have h = (7r/2 + s)/(2v^) = 1/2 as the additional waste per unit 
length of left-hand or right-hand sides. 

Now we should take a/b = 2 — \/3 as the ratio of height over width that yields the optimal 
balance between the waste along the sides of the enclosing rectangle and hence the minimum 
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area. Figure 16.11 includes a horizontal line at the level of the aspect ratio 2 — y/S. 
7. Concluding remarks 

As pointed out in the beginning of the paper, we hope that these computational results will 
lead to proofs of optimality for larger (or even infinite classes of) n. 

Our experiments suggest that the frequency of occurrence of non-hexagonal best packings 
increases with n. It is not clear whether or not the frequency has a limit as n oo and if 
it has, whether of not this limit is smaller than 1. While it is not easy to find small n > 14 
for which the best packing is not perfectly hexagonal, it might be more difficult to find large 
n for which the best packing is perfectly hexagonal. Do there exist infinitely many n, for 
which the densest packing of n circles in a rectangle is hexagonal? 

A related phenomenon is conjectured to hold for n = N{k) of the form N{k) = |(afc + 
l)(6fc + 1) where ak and bk are given by: 

ai = 1, a2 = 3, ak+2 = 4afc+i — ak 

hi = 1, 62 = 5, hk+2 = 46fc+i - h 

so that, A^(2) = 12,iV(3) = 120,A^(4) = 1512, etc. The fractions are actually (alternate) 
convergents to and it has been conjectured by Nurmela et al. |NORj that for these n, a 
"nearly" hexagonal packing of n circles in a square they describe is in fact optimal. 

In our case, one should seek alternate convergents to \/3 + 3/2 which yields sequences 
and bk given by: 

ttk = 2vk+i - Vk, bk = 2vk, 

where 

fo = 0, vi = l, Vk+2 = 4:Vk+i - Vk, k = 2,3,... 

Thus, 

(ai, 61) = (7, 2), (a2, b^) = (26, 8), (ag, 63) = (97, 30).... 

so that A^(2) = 208, A^(3) = 2910, etc. No such N{k) for = 2, 3, ... is indeed found to be 
irregular in our experiments. The best packing found experimentally for such an N{k) has 
h = bk alternating rows of full length w = ak with /i_ = 0, s = 0, in the notation of Section HI 
and Tables mU and W3 
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